NASA Technical Memorandum 107467 
AIAA-97-0150 


Collective Interaction of a Compressible 
Periodic Parallel Jet Flow 


Jeffrey Hilton Miles 
Lewis Research Center 
Cleveland, Ohio 


Prepared for the 

35th Aerospace Sciences Meeting and Exhibit 

sponsored by the American Institute of Aeronautics and Astronautics 

Reno, Nevada, January 6-10, 1997 



National Aeronautics and 
Space Administration 



Collective Interaction 
of a Compressible 
Periodic Parallel Jet Flow 

Jeffrey Hilton Miles* 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Oil 44135 


Abstract 

A linear instability model for multiple spatially peri- 
odic supersonic rectangular jets is solved using Floquet- 
Bloch theory. The disturbance environment is investi- 
gated using a two dimensional perturbation of a mean 
flow. For all cases large temporal growth rates are found. 
This work is motivated by an increase in mixing found 
in experimental measurements of spatially periodic su- 
personic rectangular jets with phase-locked screech. The 
results obtained in this paper suggests that phase-locked 
screech or edge tones may produce correlated spatially 
periodic jet flow downstream of the nozzles which creates 
a large span wise multi-nozzle region where a disturbance 
can propagate. The large temporal growth rates for ed- 
dies obtained by model calculation herein are related to 
the increased mixing since eddies are the primary mech- 
anism that transfer energy from the mean flow to the 
large turbulent structures. Calculations of growth rates 
are presented for a range of Mach numbers and nozzle 
spacings corresponding to experimental test conditions 
where screech synchronized phase locking was observed. 
The model may be of significant scientific and engineering 
value in the quest to understand and construct supersonic 
mixer-ejector nozzles which provide increased mixing and 
reduced noise. 

I. Introduction 

Interest in proving the economic and environmental 
feasibility of a high-speed civil transport has stimulated 
studies of mixing enhancement in lobed mixer-ejector 
nozzles. By enhancing mixing the ejector length can be 
reduced with the same amount of noise suppression. In 
order to obtain information on such flows simpler configu- 
rations are studied. In particular, a simple mixer nozzle 
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configuration consisting of multiple rectangular nozzles 
with a synchronized screech instability was studied by 
Taghavi and Raman 1 and Raman and Taghavi 2,3 . This 
nozzle showed increased mixing with the jets synchro- 
nized. The same behavior is shown in a study of the 
efTect of edge tones on multiple jet mixing of high-speed 
flows by Krothapalli et. al. 4 using the nozzle described 
by Krothapalli et. al. 5 . 

A study of an array of subsonic jets imbedded in a 
square network by Villermaux and Hopfinger 6 and Viller- 
niaux, Gagne, and Hopfinger 7 found the existence of 
propagating waves along the lattice. In addition, it was 
found that two jets separated by a distance up to five 
meshes had correlated oscillations. 

In this paper, the temporal dynamics produced by the 
collective interaction of jets is discussed. It is proposed 
that at some point before the jets merge it is possible to 
investigate the flow dynamics using a model based on a 
two dimensional perturbation of a mean flow. This paper 
shows a large-scale propagation of instabilities w r ith high 
growth rates may occur. 

For single nozzles a reduction in mixing and growth 
rates with increasing Mach number has been demon- 
strated experimentally by many investigators 8-12 . Corre- 
sponding linear stability analysis of single nozzles shows 
results that are similar to the experimental studies 13-15 . 
This is attributed to the fact that eddies are the primary 
mechanism that transfer energy from the mean flow to 
the large turbulent structures. However, the following 
study is based on the idea that these experimental and 
theoretical results do not apply to the mixing of multi- 
ple supersonic rectangular jets with phase locked screech. 
This paper is based on a linear stability analysis of com- 
pressible periodic parallel jet flows which was undertaken 
to obtain results related to lobed mixer nozzles. In this 
study, the lobed nozzle design concept is extrapolated in 
a one dimensional manner to arrive at an array of parallel 
rectangular nozzles separated by a distance s where the 
smaller dimension of each nozzle is xv /v and the longer di- 
mension b is taken to be infinite. Note that it is assumed 
that even widely spaced rectangular jets w hich are phase- 
locked by screech are coherent spatially at some distance 
from the nozzle. Consequently, in this linear stability 
analysis of pertubations about the mean flow, it is the 
collective behavior of compressible periodic parallel jet 
flow that determines the nozzle interaction. 

For each operating condition, the unstable wave is as- 
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sunied to grow at the maximum rate possible. In this 
paper, the behavior of the maximum growth rate solu- 
tions is discussed and only the trace of solutions for the 
group that produced the maximum growth rate for each 
case studied is presented. 


II. Results 

The nozzle configuration is shown in Figure 1. In this 
paper, the flow is compressible and the velocity profile 
perpendicular to the flow is adapted from an equation 
used by Monkewitz 16 in a study of the absolute and con- 
vective instability of two-dimensional wakes. A discus- 
sion of the problem formulation is given in Appendix A. 
A typical velocity profile is shown in Figure 2. 

The linear stability analysis is done using Floquet- 
Bloch theory. It is assumed that in the region of in- 
terest a coherent wave can propagate and that this wave 
can be described in terms of a mean flow* perturbation. 
This type of analysis has been applied by Beaumont 17 
to an incompressible flow with a sinusoidal velocity pro- 
file perpendicular to the flow. This analysis procedure is 
discussed in Appendix B. 

Stability information is obtained using the flow model 
described in Appendix A and the Floquet- Bloch method 
described in Appendix B. The flow disturbance is charac- 
terized by a real wave number, k } and a complex relative 
phase velocity , c = c r + i c*. For a given value of jet 
Mach number, Mj t ratio of inter jet spacing to rectangu- 
lar nozzle smallest dimension (s/w^), and c*, a range of k 
values are studied to determine if a growing disturbance 
characterized by a periodicity paramater T, and a con- 
vective phase velocity c P exists. The computer program 
evaluates solutions at one hundred fixed values of c r in 
the range — 1 < c r < 1 . A solutions at a given value c r 
is accepted if the calculated value of T r is smaller than 
l.E — 06 and the calculated value of 6 is less than 2 ( 
6 is defined in Appendix B). All solutions are tabulated 
and a further search is made in the c r region where T r is 
smallest to find the desired result. The reported results 
at each value of k are limited to three: no solution, one 
solution or two solutions. It is possible that more than 
two solutions exist. The model was developed to study 
wave growth, £,* = 0.5fcc,* over a range of Mach numbers 
and flow geometries for compressible periodic parallel jet 
flow when the flow is correlated between the jets. 

The Mach number and spacing for the conditions stud- 
ied correspond to cases where phase locking was achieved 
using synchronized screech by Taghavi and Raman 1 and 
Raman and Taghavi 2,3 . 

For each condition studied, solutions for a range of k 
values at a given value of c, were produced to find the 
region where the growth rate maximum, w, = O.SJtc,*, of 
the unstable wave occurred. The value of c, used were 
between 0. and 1. using steps of 0.1. The value of k 
used started at 0 and increased by 0.005. In general, 


blocks of 50 k points were examined at a one time and 
the calculation for a particular value of c, was abandoned 
if the current block of 50 points and the previous block 
of 50 points had no solutions. 

To provide some information on the solution space, the 
trace of c r and T; for the group of solutions at the value of 
c, that produced the maximum growth rate for each case 
studied is presented. Figures 3 through 17 show typical 
stability plots of phase speed c r , and T; as a function of 
wave growth, w* = O.Sifcc,*. 

The following features of these plots are noteworthy. 
Examination of the plots shows that at low Mach num- 
bers the solutions are double valued at the lower growth 
rates (Figs. 3, 4, 15 and 16). The trace of the plots of c P 
and T,- have breaks where no solution was found ( Figs. 
3, 4, 8, 13 and 15). The same Mach number (1.4) was 
used to calculate the results shown in Figure 6 ( ^=5.16 
) and Figure 17 ( ~=5 ). The slight difference in spac- 
ing produces small differences in the plots of c r and T, 
indicating that the calculations are internally consistent. 

However, the most important feature is that for the 
solution with the maximum growth rate the value of r,* 
tends to be zero. This means that the periodicity of 
the fastest growing instability wave is the same as the 
periodicity of the nozzle geometry. Consequently, it is 
possible to achieve large growth rates without an infinite 
array of nozzles. This analysis indicates an array of four 
or five nozzles should behave like an array of forty or fifty 
nozzles. 

The stability model is not related to screech. However, 
it does depend on the presence of a large span wise multi- 
nozzle region where a coherent wave can propagate. In 
this paper, it is suggested that this region can be created 
by phase locked screech or edge- tones. 

For each case, the solution values at (w,) mar are given 
in Table I which shows the fully expanded Mach number, 
A/ ; , the flow profile parameter, A, and the ratio of the 
ratio of nozzle spacing, s, to narrow width of rectangular 
nozzle, ww, ,the real and imaginary values of the phase 
speed c, the parameter T,*, the wave number it, the growth 
rate, (u>) mar , and the scale factor L *. The nozzle width, 
is 0.0069m. 

The growth rates are large. In addition, examination of 
table I indicates a tendency for the value of c, to decrease 
with Mach number. However, the corresponding value of 
k tends to increase with Mach number. Consequently, 
the maximum value of the growth rate, (w,-) mflri does 
not decrease drastically with Mach number. These large 
growth rates of spatially coherent waves might explain 
the increased mixing observed when the flow from lin- 
ear arrays of rectangular nozzles is synchronized by edge 
tones as observed by Krothapalli 4 or by screech Taghavi 
and Raman 1 , and Raman and Taghavi 2,3 since eddies 
are the primary mechanism that transfer energy from 
the mean flow to the large turbulent structures. 
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III. Concluding Remarks 


A linear instability model for a large span wise multi- 
nozzle region where the disturbance environment can be 
investigated using a two dimensional perturbation of the 
mean flow has been presented. The results indicate that 
an instability wave can occur and that this type of disur- 
brance has a large growth rate. In all the cases studied 
the most unstable wave has the same periodicity as the 
nozzle array. 

It is conjectured that multiple supersonic rectangular 
jets phase-locked by screech or subsonic jets phase-locked 
by edge tones may exhibit a high growth rates down- 
stream of the nozzles. Consequently, the model may 
explain the increase in mixing observed in multiple jets 
phase locked by screech or edge tones. 

This work was conducted with the expectation that 
multi-jets with synchronized screech could provide in- 
creased mixing and reduced aerodynmaic acoustic noise. 
The model may be of significant scientific and engineering 
value in the quest to understand and construct supersonic 
mixer-ejector nozzles 


Table I. Calculation results at (w,) mar 


M } 

A 

“■'V 

c, 

C r 

r, 

k 

( w « )m« r 

L* 

0.5 

1.24295 

5.0 

0.9 

-2 878747E-1 

97755581 E-3 

1.205 

0.54225 

3.66Q5E-4 

0.8 

1.24295 

5.0 

0.8 

-2.551930E-1 

2.955246E-2 

1.345 

0.538 

3.6605E-4 

1.25 

1.294735 

4.0 

0.6 

-2.475852E-1 

1.028344E-1 

2.06 

0.618 

5.4908E-3 

1.3 

1.28551 

4.15 

0.6 

-2.584690E-1 

3.495140E-2 

1.935 

0.5805 

5 6555E-3 

1.35 

1.22414 

5 5 

0.6 

-1.648727E-1 

3.762793E-2 

2.06 

0.4785 

7.1380E-3 

1.4 

1.24295 

5.0 

0.6 

-2 084984E-1 

7.771379E-2 

1.41 

0.423 

3.6605E-4 

1.4 

1.23649 

5 16 

0.6 

-1 .860186E-1 

1.184881E-1 

1.935 

0.411 

6.7647E-3 

1.45 

1.173897 

7 5 

0.5 

-2 394745E-2 

2.976108E-2 

2.060 

0.515 

9 3344E-3 

1.5 

1.19688 

6 43 

0.5 

-9 262075E-2 

6. 455086 E-2 

1.91 

0.411 

8.1593E-3 

1.55 

1.139382 

10 0 

0.4 

8.28761 2E-2 

9. 056775 E-2 

2 245 

0.449 

1.2079E-2 

1.6 

1.160343 

8.32 

0.4 

9 .90190 1 E-3 

1 . 042738 E-l 

2.240 

0.448 

1.0234E-3 

1.65 

1.1258595 

11.5 

0.3 

8.145849E-2 

3.5661 17E-2 

2.55 

0.3825 

1.3727E-2 

1.7 

1.14881 

9 16977 

0.3 

1 .51041 1 E-2 

9.125563E-2 

2.615 

0.39225 

1.1168E-2 

1.75 

1.111395 

13 7 

0 3 

1 .772602E- 1 

9. 963684 E-3 

1 89 

0.2835 

1 .6 143E-2 

1.8 

1.13665 

10.27 

. 0.3 

7.961790E-2 

4.183396E-3 

2.095 

0.31425 

1.2376E-2 
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Appendix A: Formulation of the problem 


_ 

m — 


mf 


Let (t/(y),0,0) be the velocity of a steady plane- 
parallel flow, where the x-axis is in the direction of the 
flow and 

U(v) = U + ^-h{y) 

where U\ is the velocity outside the jet, U 2 is the mean 
centerline jet velocity, U = A U = U 2 - U\ t and 

h(y) is the velocity profile function which varies from -1 
to 1. 

The flow field is perturbed by introducing wave distur- 
bances in the velocity and pressure with amplitudes that 
are a function of y. Thus, 

(u,v, w,p) 

= («(y)> *>(y). *(y), p(y)) exp |; (ix + iz -urj J . 

Where 


k = kL* , 
i = CL * , 

. _ uL*_ 
U ~ A U’ 

U) _ u> 

I - T 


kAU 

and we define c as follows 


AU 


c = 


AU 


U c 
AU + 2' 


By definition k is real positive number that represents 
the wavenumber in the z-direction, i is the wavenumber 
in the ^-direction, c r is the relative phase velocity, and 
^ is the amplification rate of the disturbance. 

From the equations of motion if nonlinear and viscous 
terms are neglected one can obtain an equation for the 
y-component of the perturbation velocity as follows: 

T' A ' 

v n - 1 >'( 4 - + 

v T A 1 


h" t T A\ h' 

(h-c) K t A } {h-c) 


= 0 (Al) 


where the primes denote differentiation with respect to 


y 

p 

A — —ik — t-r- + m 2 ik 

k 


and from Crocco’s Equation 18 

t(y) = 

r 3 (1 + %)) r 2 

Ti 2 K Ti ’ 


“o( m i) 2 (7- 1) 


(%) + !)(%)-!) 

4 


where 


mi = 


AU A U y/Ti 


s/% 




a 2 V% m *v*T 


In this paper, the velocity profile function, h(y), is pe- 
riodic such that 


h(y + 2 tt) = h(y). 

The velocity profile h(y) is not any exact solution of 
the Navier-Stokes equation, but it can be considered as 
a simple model of some real periodic flow. 

The velocity profile h(y) discussed herein is given by 

My) = 1 - 2/(y) 

where the function /(y) is given by 

1 


/(v) = 


1 + ( sinh (^7Tj)) 


18 


7 / — A(— 1 H — ), 

7 T 

values of A for the ratios of ( s/w //) used herein are 
given in Table II and y goes from 0.0 to 2 tt. The 
profile function f(y) is adapted from an equation used 
by Monkewitz 16 in a study of the absolute and con- 
vective instability of two-dimensional wakes. Only two- 
dimensional disturbances will be considered. A schematic 
of the nozzle geometry is shown in Figure 1. A typical 
velocity profile using A = 1.5 is shown in Figure 2. 


A' = 2m 2 ik—~ 
4 
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Appendix B: Floquet-Bloch theory 


Since 


Since the basic flow velocity profile , /(y), is peri- 
odic, equation (Al) is an example of a Floquet-Bloch 
problem. The mathematics of solving Floquet-Bloch 
type problems is discussed by Ince 19 , Hochstadt 20 , and 
Zwillinger 21 . Applications to solid state physics are dis- 
cussed by Sachs 22 , Brillouin 23 , and Dekker 24 . Applica- 
tions to spatially periodic flow are discussed by Lorenz 25 , 
Green 26 , Beaumont 17 , and Gotoh 27,28 . 

The paper by Beaumont 17 and the description of the 
Floquet-Bloch theorem by Hochstadt 20 were particularly 
useful in guiding this research. 

A survey of the spatially periodic flow literature is pre- 
sented by K. Gotoh and M.Y. Yamada 29 . 

The second order differential equation can be described 
by a system of first order differential equations. Let 

v — Xi 
v = 2*2 


so that Eq. Al can be rewritten as the system 


where 


and 



T' A* 
T A 


(Bl) 


D -- 


Tj rr 4* Aik — (— H — — ) — — 

(h-c) v T A (h — c) 


If <$(y) is a fundamental matrix solution of equation ( 
(Bl ) ) such that 


$(0) = I 


where I is the identity matrix, then from the Floquet- 
Bloch theorem 


$(y + 2tt) - $(y)<I>(27r) 


We now introduce two solutions of equation ( Bl ) 
with initial values at y = 0.0 


m = 


<M0) ^(0) 

>i'(0) ^'(0) 



'l o' 


0 1 


Next we seek the eigenvalues of $( 27 t) 

l*(2*) - ^I| 

_ <£i( 2 jr)-/i <j> 2(2t) 

- <j>\'{ 2 v) 4>2{2ir) - /i 

= /i 2 - (4> i(27t) + 2n))fi 
+(^i(27r)^ 2 '(27r) - 
^(2^/(2^)) 

= 2jt) + <A 2 '(27r))/r + 1 = 0 (B2) 


^i(2jt)^ 2 , (2jt) - ^ 2 (27r)^ 1 '(2w) = |<K(2r)| = |<*(0)| = 1 

The independent solutions of equation (Bl) have the 
form 


<t> = *(y)exp(^j^y) = X(y)exp(ry) 

The parameter T specifies the period of the eigenfunc- 
tion 4>. If T is real the eigenfunction grows or decays at 
infinity. Consequently, only imaginary values of F are 
acceptable. Thus the eigenfunction oscillates in space 
and is called a continuous mode. The disturbance with 
r t = 1/n, where n is a nonzero integer, has a period 2nn. 
One with T, = 0 has the same period 27 t as the main flow, 
while an irrational value of F,* means the disturbance is 
aperiodic. Note that the parameter F does not appear 
in the flow equation, but is due to the Floquet-Bloch 
theorem. 

Solutions of Bl are thus of the form 

X l (y + 2x) = n l X l (y) 

V 2 (y + 2jr) = y .2 X 2 (y) 

where hi and fi 2 represent the zeros of (B2), provided 
they are distinct. 

In general, these solutions will not be periodic. 

Conditions for periodic solutions can be found as fol- 
lows 

Let pi = e 1 * 1 and H2 = 

Then from equation (B2) 

cos(^i) = <2 i(2tt) -f <^> 2 ^( 2 ^) = 6/2 

Consequently, for a solution to be periodic 6 must be 
real and \6\ smaller than 2. 

The constants /i are termed the characteristic multi- 
pliers of the Floquet-Bloch system (Bl) and the corre- 
sponding characteristic exponents are determined by the 
relation T = T r + «T,- = lo ^ = |jr + ifp. 
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range of Mach numbers and nozzle spacings corresponding to experimental test conditions where screech synchronized 
phase locking was observed. The model may be of significant scientific and engineering value in the quest to understand 
and construct supersonic mixer-ejector nozzles which provide increased mixing and reduced noise. 
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